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Abstract 

We examine the effects of a periodically varying flow velocity on the standing and travelling 
wave patterns formed by the flow-distributed oscillation (FDO) mechanism. In the kinematic 
(or diffusionless) limit, the phase fronts undergo a simple, spatiotemporally periodic longitudinal 
displacement. On the other hand, when the diffusion is significant, periodic modulation of the 
velocity can disrupt the wave pattern, giving rise in the downstream region to travelling waves 
whose frequency is a rational multiple of the velocity perturbation frequency. We observe frequency 
locking at ratios of 1:1, 2:1 and 3:1, depending on the amplitude and frequency of the velocity 
modulation. This phenomenon can be viewed as a novel, rather subtle type of resonant forcing. 
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Interest has been mounting recently on mechanisms of pattern formation in open reactive 
flows. The combination of reaction, advection and diffusion, together with the effect of 
an upstream boundary condition, leads to mechanisms such as flow-distributed oscillations 
(FDO) Q-j^, a general category of stationary patterns referred to as "flow and diffusion- 
distributed structures" (FDS) Q- 3> and the differential flow instability (DIFI) Q 11 1. 

Our focus here is on FDO. Due to the equivalence of flow in reaction- 

advection-diffusion (RAD) systems and linear growth of the spatial domain of a reaction- 
diffusion system, and the existence of cellular oscillations in segmenting tissue, FDO was 
shown to be involved in the axial segmentation occurring^during biological development 



18 1 . Given the pulsating growth of certain organisms |2l|-|22|, including human embryos 
23[, we study here the consequences of a periodically modulated flow v{t) on FDO. 

The systems of interest are described by the RAD equation without differential transport: 



where D is the diffusion constant, v(t) is the flow velocity, U(x, t) is an iV-dimensional vector 
of dynamical variables (concentrations of species) and the local dynamics given by the vector 
valued rate function f(U) has an attracting limit cycle. If v (t) is constant, this system can 
support flow-distributed oscillations (FDO) controlled by the upstream boundary condition 
U(0, t). In the simplest constant boundary condition sets the phase of each oscillat- 

ing fluid element as it enters the medium, and the periodic recurrence of the same phase as 
the fluid travels downstream results in stationary waves. Oscillating boundary conditions 



result in travelling waves. 



id ] Diffusion can modify the effective dynamics of the 



medium as it travels downstream and even extinguish the oscillations. |2J [2| Equation (1) 
is also relevant to media such as linearly growing organisms, as it can be reinterpreted by 
means of a Galilean transformation as rep resenting a stationary medium with a boundary 
(the growth tip) moving at speed v (t) . [3 [3 13 13 

We examine the effect of a sinusoidally varying velocity v(t) — Vq + 5 v cosu v t. In the 
kinematic limit of vanishing diffusion D/v 2 — > the wave pattern undergoes a simple, 
calculable longitudinal displacement which is periodic in both time and space. Away 
from the kinematic limit, however, we observe a type of nonlinear resonance. Relatively 
small disturbances of an FDO wave pattern are magnified with downstream distance until 
the wavefronts break. This rupture generates travelling waves in the downstream region 
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whose temporal frequency is a rational multiple of the velocity perturbation frequency. We 
observe 1:1, 2:1 and 3:1 ratios depending on the frequency and amplitude of the velocity 
perturbation. 



For the numerical examples, we use the FitzHugh-Nagumo-type (FN) dynamics 24] 



e(X -X 3 -Y)\ 
f(X,Y)= (2) 
y-Y + aX + P J 

with e = 5, a = 2 and (3 = for the local rate function. At these parameter values, the 
local system has a limit cycle and a moderately strong nonlinearity. The frequency of the 
limit cycle oscillation is u « 27r(0.43). In all simulations, we set D = 1 and vary only v 
and S v . 

When the ratio D/v 2 is sufficiently small, diffusion is relatively unimportant and each 
individual fluid element behaves approximately as an independent oscillator obeying the 
local dynamics whose initial phase is set by the boundary condition as it enters the flow 
from the upstream end. j^] The phase fronts can then be calculated by pure kinematics: the 
oscillation phase of a fluid element at a particular time and location depends on its initial 
phase when it entered the flow and how long ago it entered the flow. For the case of a 
stationary boundary condition (i.e., constant phase at the boundary) the result is that the 
location of the phase front for a particular value of the oscillation phase <j) is given by 

x[<p,t) = Vq 1 {smu v t — sm.u v (t )), (3) 

UJq U) v U>q 

where ujq is the frequency of the local oscillator. When S v = 0, this reduces to the 
simple linear mapping between position and phase that characterizes stationary FDO waves. 
Note that: 1) The amplitude of the displacement of the phase fronts (second, time- 
dependent term in eq. EJ) depends on 5 v /u v , implying that faster (higher frequency) 
velocity modulation has less of an effect than slower modulation. 2) The displacement 
is periodic in time with the same frequency as the velocity perturbation, but 3) it is also 
periodic in and thus in space. The fronts for which is a multiple of huoq/uj v are not 
displaced, and they occur periodically at positions nnv /u v . The spatial periodicity can be 
understood by considering the trajectories of fluid elements entering the system at different 
times. Different elements enter at different points in the velocity modulation period and 
thus begin their downstream travel at different initial velocities. Over any multiple of the 
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FIG. 1: Periodic modulation of standing and travelling FDO waves in the kinematic limit. Note 
the periodicity in both time and space. In each space-time diagram, the trajectory of a co-moving 
point is shown as a guide to the eye (dashed white line). (A) Stationary waves: vq = 8, S v = 2, 
u v = 27r(0.1) ~ wo/4. Some wavefronts remain stationary while others wiggle back and forth. (B) 
Upstream travelling waves with wave frequency u>tw = 27r(0.25) subject to a modulated velocity 
field with v = 5, 5 V = 2, uj v = 2vr(0.02). 



modulation period, however, the velocity averages to v and thus all elements reach the 
same position at the same phase when one full period has passed, regardless of when they 
started. 

When the boundary condition is oscillatory instead of stationary, travelling waves are 
generated. Just as in the stationary case, velocity modulation causes a periodic longitudinal 
displacement of the travelling wavefronts. The modulation of stationary and travelling 
waves is illustrated in figure ??. 

When diffusion is unimportant, neighboring fluid elements do not interact and the be- 
havior of FDO patterns can be explained by pure kinematics. Each co-moving fluid element 
follows the limit cycle defined by the batch reactor dynamics. However, significant diffusion 
alters the dynamics. The flow velocity modulation then introduces a periodic variation in 
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the local environment of each fluid element. The strength of the local gradient is differ- 
ent for each co-moving element and diffusion therefore affects the dynamics differently at 
different locations. This differential effect can magnify the small kinematic effect of the 
velocity perturbation, leading to larger differences in the dynamical variables and eventually 
to disruption of the smooth FDO waveforms. 

Figure El shows several examples of this phenomenon, in which quite subtle modulations 
of a stationary wave pattern become magnified with increasing downstream distance and lead 
to the breaking and reconnection of wave fronts in the downstream region. The simulations 
in these examples were all done at an average flow velocity of Vq = 3. By comparison, the 
boundary between absolute and convective instability of the Hopf/FDO instability occurs at 
vac ~ 2.82. Stationary waves controlled by the boundary remain possible at velocities well 
below this threshold, however. Thus, while the flow velocity is far from the kinematic 
limit, it lies well within the regime where boundary-controlled stationary waves are stable 
in the absence of flow modulation. The minimum velocity Vq — 5 V never falls below vac 
except in figure EK, and then only by a small amount. 

In figure EK, the perturbation frequency is very close to the natural frequency of the 
chemical oscillator. The effect of the flow modulation is visible in this space-time plot as 
a slight pulsation of the wavefronts. The pulsation becomes stronger at positions farther 
downstream, until there is a transition to a region of nearly uniform synchronous oscillation, 
synchronized to the period of the flow modulation. If the modulation frequency is changed, 
the pattern in the downstream region remains synchronized to the modulation, and the result 
is either upstream or downstream travelling waves. An example of upstream waves is shown 
in figure |2J3- In this case the velocity modulation is at a frequency lower than the intrinsic 
natural frequency of the medium. As one can see from the figure, the system's response 
to the velocity perturbation is nonlinear. Instead of a simple sinusoidal displacement, the 
stationary wavefronts develop a series of sharp cusps. At a certain downstream position, 
the wavefronts break and reconnect, and the periodic disturbances become the source of 
a set of travelling waves with frequency equal to the modulation frequency, just as if the 
boundary were being driven at that frequency. 

Just as in the case of ordinary FDO phase waves driven by a perturbation at the boundary 
DQD perturbations slower than the intrinsic frequency give rise to upstream travelling 
waves. (In general, the phase velocity, wavelength and frequency of the travelling waves 
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FIG. 2: Examples of the breakup of waves due to a velocity perturbation. All examples have 
average flow velocity vo = 3. The dotted white line in each frame represents vq. A) uj v = 
27r(0.43) ps loq, S v = 0.215. A slight pulsation with frequency u v is visible in the standing 
wavefronts, baecoming more pronounced downstream. Near x = 100 there is a transition to 
uniform oscillation. B) u v = 27r(0.2) ps ujq/2, S v = 0.1. Periodic disturbance of the standing 
waves becomes sharper with increasing distance, and there is a transition to travelling waves with a 
1:1 frequency ratio. Near the transition these travelling waves propagate in a saltatory manner but 
they grow smoother with further downstream distance. C) u v = 27r(0.15) ps wo/3, 5 v = 0.1125. 
As in (B) there is a transition to travelling waves near x = 100, but in this case the waves do 
not smooth out with downstream distance. Instead, there is a second transition at x=200 to 
waves with twice the velocity modulation frequency (2:1 resonance). D) u v = 27r(0.10) ps wo/4, 
5 V = 0.05. A series of transitions leads to waves with three times the perturbation frequency (3:1 
resonance) . 6 



obey the kinematic relationships discussed in 

and Q .) These waves propagate 
irregularly in the region just downstream from the transition, but with increasing down- 
stream distance they become smoother. The temporal frequency of the waves is locked 
to the velocity modulation frequency. Figure |2P shows a more complicated situation 
with two consecutive transitions. The first transition to travelling waves occurs much as in 
133. However, instead of smoothing out with downstream distance, these waves propagate 
irregularly and develop a second instability at a position farther downstream, leading to 
travelling waves with a temporal frequency exactly twice that of the velocity perturbation. 
This can be viewed as a form of 2 : 1 frequency locking. The latter travelling wave smooths 
out with downstream distance and appears to be the final asymptotic waveform. An 
asymptotic waveform at three times the velocity perturbation frequency is also possible, as 
in figure Hp. In general, a sequence of transitions leads to successive regions of station- 
ary, 1:1, 2:1, 3:1, etc. waves. In the particular case of fig. EP, the three transitions are 
quite close together. Which asymptotic waveform is selected, and the exact distances from 
one transition to the next, depend in nontrivial ways on the perturbation frequency and 
amplitude. We will explore this dependence in subsequent work; a behavior somewhat 
analogous to Arnold tongues seems to occur. We have observed asymptotic waves in 1:1, 
2:1 and 3:1 ratios to the perturbation frequency, but we have not yet observed other rational 
multiples such as 2/3. Interestingly, the tendency of stationary waves to break is strongest 
not at the intrinsic natural frequency of the oscillator, but at approximately 0.65o;o- 

Immediately downstream from any transition point, the travelling waves generally prop- 
agate with a pulsating phase velocity, but become smoother with increasing distance down- 
stream. Such behavior was observed in both experiments and numerical simulations for 
waves forced at the boundary under a steady flow velocity [4J. In that case, the pulsating 
phase velocity was due to a mismatch between the oscillations driving the waves and the 
limit cycle of the intrinsic dynamics in the flow reactor. The explanation in this case is 
the same. Instead of being driven by an oscillation at the inflow boundary, however, these 
travelling waves are driven by an oscillation induced by the flow velocity modulation. As 
in the case without flow modulation, diffusion tends to smooth the jumping waves as they 
travel downstream, unless the velocity perturbation induces a second instability as in figures 
I2P,D. 

The spatiotemporal resonance manifested in wave-front disruption and frequency locking 
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is novel. While previous studies of spatiotemporal resonance involved direct, global per- 



turbations of the local dynamics 26j|27l|. in the present case the perturbation acts only at 
the inflow boundary. This becomes more evident when one considers the equivalent grow- 
ing reaction-diffusion system in the co-moving frame ^|^|, where the velocity does not 
enter into the dynamical equations except via the boundary condition. Yet this boundary 
effect propagates into the spatial domain where it leads the breakup of waves and resonant 
frequency locking. 



Due to the above-mentioned equivalence of flow and growth 
nomenon should be observable in experiments such as those of Q| 



I.17L18] 



the same phe- 
which use a station- 



ary medium with a moving boundary, if the velocity of the boundary is modulated. It may 

in which growth is pulsatile. 



also be relevant to biological situations |21|- 



[1] S.P. Kuznetsov, E. Mosekilde, G. Dewel, and P. Borckmanns, J Chem. Phys. 106, 7609 (1997). 
[2] P. Andresen, M. Bache, E. Mosekilde, G. Dewel and P. Borckmanns, Phys. Rev. E. 60, 297 

(1999) . 

[3] M. Kaern and M. Menzinger. Phys. Rev. E. 60, R3471 (1999); 62, 2994 (2000); P. Andresen, 

E. Mosekilde, G. Dewel and P. Borckmans, Phys. Rev. E 62, 2992 (2000). 
[4] M. Kaern, M. Menzinger. Phys. Rev. E. 61, 3334 (2000). 
[5] M. Kaern and M. Menzinger M., J. Phys. Chem. 106, 4897 (2002). 

[6] J.R. Bamforth, S. Kalliadasis, J. H. Merkin, S. K. Scott, Phys. Chem. Chem. Phys. 2, 4013 

(2000) . 

[7] J.R. Bamforth, J. H. Merkin, S.K. Scott., R.Toth and V. Gaspar, Phys. Chem. Chem. Phys. 
3, 1435 (2001). 

[8] J.R. Bamforth, R. Toth, V. Gaspar and S.K. Scott, Phys. Chem. Chem. Phys. 4, 1299 (2002). 

[9] P.N. McGraw and M. Menzinger, Phys.Rev. E 68 066122 (2003). 
[10] A.B. Rovinsky and M. Menzinger, Phys.Rev.Let. 69, 1193 (1992), 70, 778 (1993). 
[11] M. Menzinger and A.B. Rovinsky, "The Differential Flow Instability," in R. Kapral and K. 

Showalter, Chemical waves and patterns. (Kluwer, Dorderecht,1995). 
[12] R.A. Satnoianu and M. Menzinger, Phys. Rev. E. 62, 113 (2000). 
[13] R. Satnoianu, P.K. Maini and M. Menzinger., Physica D 160, 79, (2001). 



S 



[14] R.A. Satnoianu, Phys. Rev. E 68, 032101 (2003). 
[15] P.N. McGraw and M. Menzinger, preprint |nTm.PS /04100l"7| (2004). 
[16] M. Kaern, M. Menzinger, A. Hunding, Biophys. Chem. 87, 121 (2000). 
[17] M. Kaern, M. Menzinger, A. Hunding, J. Theor Biol. 207, 473 (2000). 
[18] M. Kaern, M. Menzinger, R. Satnoianu and A. Hunding, Faraday Disc. 120, 295 (2002). 
[19] M. Kaern, R. Satnoianu, A. Munuzuri and M. Menzinger, Phys. Chem. Chem. Phys. 4, 1315 
(2002). 

[20] M. Kaern, D.G. Miguez, A. P. Munuzuri and M. Menzinger, Biophys. Chem. 110, 231 (2004). 
[21] R. Lopez Franco, S. Bartnicki Garcia, C.E. Bracker Proc.Natl.Acad.Sci.USA 91, 12228 (1994). 

M. Messerli, K.R. Robinson. Journal of Cell Science 110, 1269 (1997). M. Messerli, K.R. 

Robinson. Plant Journal 16, 87 (1998). M. Messerli, G. Danuser, K.P. Robinson. Journal of 

Cell Science 112, 1497 (1999). 
[22] D.A. Care, J.R. Crush. New Zealand Journal of Agricultural Research 47, 67 (2004). 
[23] M. Lampl, P. Jeanty. American Journal of Human Biology 15, 667 (2003). 
[24] C.B. Muratov and V.V. Osipov, Phys. Rev. E 54, 4860 (1996). S.P. Dawson, M.V. D'Angelo 

and J.E. Pearson, Phys. Lett. A 265, 346 (2000). A. Hagberg and E. Meron, Chaos 4, 477 

(1994). 

[25] PV. Kuptsov, Physica D 197, 174 (2004). 

[26] F. Sagues, D.G. Miguez, E.M. Nicola, A. P. Munuzuri, J. Casademunt, L. Kramer, Physica D 
199, 235 (2004). 

[27] J.H.P. Dawes, CM. Postlethwaite, M.R.E. Proctor, Physica D 191, 1 (2004). J. Porter, E. 
Knobloch, Physica D 159, 125 (2001). 



9 



